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THE SELBERG ZETA FUNCTION FOR CONVEX CO-COMPACT 

SCHOTTKY GROUPS 

LAURENT GUILLOPE, KEVIN K. LIN, AND MACIEJ ZWORSKI 


Abstract. We give a new upper bound on the Selberg zeta function for a convex co-compact 
Schottky group acting on in strips parallel to the imaginary axis the zeta function is 

bounded by exp(C|s|^) where <5 is the dimension of the limit set of the group. This bound 
is more precise than the optimal global bound exp(C|s|^'*'^), and it gives new bounds on 
the number of resonances (scattering poles) of The proof of this result is based on 

the application of holomorphic L^-techniques to the study of the determinants of the Ruelle 
transfer operators and on the quasi-self-similarity of limit sets. We also study this problem 
numerically and provide evidence that the bound may be optimal. Our motivation comes 
from molecular dynamics and we consider as the simplest model of quantum chaotic 

scattering. The proof of this result is based on the application of holomorphic L^-techniques to 
the study of the determinants of the Ruelle transfer operators and on the quasi-self-similarity 
of limit sets. 


1. Introduction 


In this paper we give an upper bound for the Selberg zeta function of a convex co-compact 
Schottky group in terms of the dimension of its limit set. This leads to a Weyl-type upper bound 
for the number of zeros of the zeta function in a strip with the number of degrees of freedom 
given by the dimension of the limit set plus one. We also report on numerical computations 
which indicate that our upper bound may be sharp, and close to a possible lower bound. 

Our motivation comes from the study of the distribution of quantum resonances - see 
for a general introduction. Since the work of Sjostrand |M| on geometric upper bounds for the 
number of resonances, it has been expected that for chaotic scattering systems the density of 
resonances near the real axis can be approxim atel y given by a power law with the power equal to 
half of the dimension of the trapped set (see ( 0 ) below). Upper bounds in geometric situations 
have been obtained in Q and 1^ . 

Recent numerical studies in the semi-classical and several convex obstacles settings, [llT|,|ll^ 
and respectively, have provided evidence that the density of resonances satisfies a lower 
bound related to the dimension of the trapped set. In complicated situations which were studied 
numerically, the dimension is a delicate concept and it may be that different notions of dimension 
have to be used for upper and lower bounds - this point has been emphasized in fl^ . 

Generally, the zeros of dynamical zeta functions are interpreted as the classical correlation 
spectrum |^ . In the case of convex co-compact hyperbolic quotients, X = r\]HI”+^ quantum 
resonances also coincide with the zeros of the zeta function - see . The notion of the dimension 
of the trapped set is also clear as it is given by 2(1 -|- S). Here S = dim A(r) is the dimension of 
the limit set of F, that is the set of accumulation points of any F-orbit in A(F) C 
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Hence we may expect that 

(1.1) ^ TOr(s) ~ , 

I Im s|<r , Re s> —C 

where mr(s) is the multiplicity of the zero of the zeta function of F at s. 

Referring for definitions of Schottky groups and zeta functions to Sections || and ^ respectively 
we have 

Theorem. Suppose that T is a convex co-compact Schottky group and that Zr{s) is its Selberg 
zeta function. Then for any Cq > 0 there exists Ci such that for | Re s| < Co 

(1.2) |Zr(s)| < Ciexp(Ci|s|'^), 5 = dim A(r). 


The proof of this result is based on the quasi-self-similarity of limit sets of convex co-compact 
Schottky groups and on the application of holomorphic L^-techniques to the study of the deter¬ 
minants of the Ruelle transfer operators. 



Figure 1. The plot of log log A^(fc)/log fc — 1, where N{k) is the number of of 
zeros with |Ims| < k, for a Schottky reflection group with S ~ 0.184. Different 
lines represent different strips, and the thick blue line gives S. 

If we use the convergence of the product representation (0 of the zeta function for Re s large 
and apply Jensen’s theorem we obtain the following 

Corollary 1. Let mr{s) be the multiplicity of a zero of Zr at s. Then, for any Co, there exists 
some constant Ci such that for r > 1 

(1.3) ^{mr(s) : r < | Im s| < r-|-1, Re s > —Cq} < Cir"*, 

where 5 = dim A(r). 
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We can apply the preceding results to Schottky manifolds: a hyperbolic manifold is said 
Schottky if its fundamental group is Schottky. The case of surfaces being of special interest: any 
convex co-compact hyperbolic surface is Schottky. With the description of the divisor of the zeta 
function through spectral data established by Patterson and Perry ||^, we can reformulate the 
preceding corollary nicely in the resonance setting. We do it only for surfaces (see below for short 
comments on higher dimensions) 


Corollary 2. Let M be a convex co-compact hyperbolic surface, Sm be set of the scattering 
resonances of the Laplace-Beltrami operator on M and mM{s) be the multiplicity of resonance 
s. Then, for any Cq, there exists some constant Ci such that for r > 1 

(1.4) ^^{ut-m(s) : s G TZ,r < \ lms\ < r -I- 1, Res > —Cq} < Cir'^, 

where 2(1 -|- 5) is the Hausdorff dimension of the recurrent set for the geodesic flow on T*M. 

This corollary is stronger than the result obtained in [ psf where the upper bound of the type 
dO) was given. In fact, the upper bound (0) is what we would obtain had we had a Weyl law 
of the form r^^^ with a remainder 0{r^). That local upper bounds of this type are expected 
despite the absence of a Weyl law has been known since [^. 

Section ^ deals with numerical computations of the density of zeros. They show that (0) 
may be true|^. In fact, in the range of Im s used in the computation we see that the number 
of zeros grows fast. If the range of Res is large (and fixed) we need very large Ims to see the 
upper bound of Corollary 2. The computations also show that our bound on the zeta function 
is optimal. For values of Zr(s) with Res negative we see that we need very large Ims to see 
the onset of the upper bound. That is not surprising since we recall in Proposition 3.2 that 
log |^r(s)| = C(|s|"+^), and that this bound is optimal (and of course 6 < n). 

We refer to Section ^ for the details and present here two pictures only. We take for P a group 
generated by compositions of reflections in three symmetrically spaced circles perpendicular to 
the unit circle, and cutting it at the angles 30° (see Fig.|| for the 110° angle). Figj^ shows the 
density of zeros of Zr in that case and Figplots the values of log | log |^r||- 


2. Schottky groups 

The hyperbolic geometry on and the conformal geometry on its boundary at infinity 

i9]HI"+^ = S” share the same automorphism group: the isometry group Isom(]HI"’+^) and the 
conformal group Conf(S") (with the conformal structure given by the standard metric on S" of 
curvature -1-1) are isomorphic. In particular any isometry g of induces on S" a conformal 

map 7 , whose conformal distorsion at the point u> G S" will be denoted by ||£) 7 (w)||. There is 
also a correspondence between balls V and spheres C on S" (for n = 2, the original setting for 
Kleinian groups, these are discs and circles) and half-spaces P and geodesic hyperplanes H in 
jjn-i-i: 25 — 'Pn9]HI"+^ and C — . Given a hyperplane Ti (a sphere C resp.), its interior 

hyperplane (ball) will be given by the choice of a component of \ Ti (§"■ \ C resp.) 

Let us review definitions of a Schottky group (see |l^, |^, and references given there), de¬ 
fined originally in conformal geometry on S'^. To the configuration of mutually disjoint geometric 
balls on the sphere S'", i = 1,... £, we associate the Schottky marked reflection group, 

r( 2 ?i,..., 2 ?£). 


^Strictly speaking, the numerical computations are done for a slighly different zeta function of conformal 
dynamical system 
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Figure 2. Density of values of log | log |Zr||/log |s| for a Schottky reflection 
group with S ~ 0.184. 


defined as the group generated by the inversions Hi in dVi^i = . The corresponding 

hyperbolic group is the Schottky marked reflection group 

for the collection, Vi,... ,Ve, of mutually disjoint hyperbolic half-spaces Vi with boundaries at 
infinity given by I?i’s. The group is generated by the hyperbolic symmetries Si in the hyperplane 
dVi (with infinite boundary dVi), i = 1,.. 

More generally, to a configuration of 2£ disjoint topological balls 

V,, i = l,...,2£, 

with a set of conformal maps "/i,i = satisfying 

7,(S"\D0=^, 

we associated the Schottky marked group 

generated by the 7 i, j = 1 ,..., £. 

We also introduce maps 

le+j =77^J = 

which map the exteriors of Ve+j ’s onto Vj. 

The Schottky group F is a free group on the £ (free) generators 71 ,..., 7 ^. If the topological 
balls T>i,i = 1,... ,2£ are geometric balls, the marked Schottky group is said to be classical. A 
classical Schottky marked group r(X>i,..., 7 i, ■ ■ ■, 7 ^), subgroup 0001 ( 8 "), can be presented 
as an isometry group of the hyperbolic space 
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Let Vi be the hyperbolic half-space with boundary at infinity given by Vi. If gi is the hyperbolic 
isometry of EI”+^ with action at infinity given by 7 i, then the classical Schottky group 

r = T[Vi,... ,X’2^,7i, ... ,7 ^), 

has a hyperbolic marking 

r = T{Vi ,.. .,V 2 i,gi, ...,gi). 

The Schottky domain \ uf^iVi is a fundamental domain for the action of T on 

A group is said to be a Schottky (reflection) group if it admits a presentation induced by 
a configuration of balls as described above above. We remark that the subgroup of positive 
isometries of a Schottky reflection group is a classical Schottky group, with a marking 

(2 ?i,... Vi-i,aiVi ,..., aiVi-i,gi,... ,gi-i ), 

where gi = aiUi. Such a group was called symmetrical by Poincare |Q. 

An oriented hyperbolic manifold M is said to be (classical) Schottky if its fundamental group 
(realized as a discrete subgroup of Isom^(IHI"+^)) admits a (classical) Schottky marking. 



Figure 3. Tesselation by the group, Tg , 9 = 110°, generated by symmetries in 
three symmetrically placed lines each cutting the unit circle in an 110 ° angle, 
with the fundamental domain of its Schottky subgroup of direct isometries, F^, 
and the associated Riemann surface The dimension of the limit set is 

S = 0.70055063.... 

A Schottky group is convex co-compact if in addition the closures Vi are mutually disjoint. 
In the convex co-compact case all non elliptic elements of F are hyperbolic. That means in 
particular that for 7 G F there exists a G Isom(]HI"+^) such that, in the Poincare model ~ 

= R+ X 

(2.1) a~'^-fa{x,y) = , (a:,y)GlR++\ 0(7) G 0(n), £(7)70. 

If P C Isom'''(]HI"+^), the conjugacy classes of hyperbolic elements, 

{71} = {72} 3 /3 G P /37i/3”^ = 72 , 

are in one-to-one correspondence with closed geodesics of r\]HI"+^. The primitive geodesics 
correspond to conjugacy classes of primitive elements of F (that is, elements which are not non¬ 
trivial powers). The magnification factor exp(^( 7 )) in & gives the length £( 7 ) of the closed 
geodesic. 

The limit set, A(r) of a discrete subgroup, P, of Isom(]HI”+^), is defined as the set in 11"'+^ = 
jjn+i y of accumulation points of any F-orbit in 11"+^: the limit set A(r) is included 

in the boundary In the convex co-compact case it has a particularly nice structure; 

furthermore, for Schottky groups, it is totally disconnected and included in = uf^iVi. The 
aspects relevant to us come from the work of Patterson and Sullivan - see and references 
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Figure 4. A typical limit set for a convex co-compact Schottky group, F C Isom(]HI^). 

given there. As will be discussed in more detail in Section^, the limit set has a quasi-self-similar 
structure and a finite Hausdorff measure at dimension S = <5(r). 

The limit set is related to the trapped set, K, of the usual scattering Q, that is the 
set of points in phase space such that the trajectory through that point does not escape to 
infinity in either direction.: if A is the diagonal oi K x K and tt* the projection from y*]H["+i on 
the trapped set K is the union of the projections Tr*{C^n) where C^n is the geodesic 
with extremities ^ and 77, both in the limit set A(r). In particular, we have 

dim AT = 2{S + 1), 

see [Q. 

To stress the connection to closed orbits let us also mention that, generalizing earlier results 
of Guillope and Lalley Perry ^ showed that 

tl{{7} : 7 primitive, £(j) < r} ~ . 

or 

3. Properties of the Selberg zeta function 
For F, a discrete subgroup of Isom(H”+^), the Selberg zeta function is defined as follows 
(3.1) Zr{s) = n n (1 “ . 

{7}“GNJ 

Here, 7 G F are hyperbolic, exp(^( 7 ) -|- iOj{'j j) are the eigenvalues of the derivative of the action 
of 7 on S" at the repelling fixed point (exp(i 0 j( 7 )) are the eigenvalues of the isometry 0 ( 7 ) in 
the normal form (^^)) and 

{7} = the conjugacy class of a primitive hyperbolic element 7. 

An element is called primitive if it is not a non-trivial power of another element. 
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We define the following map on T) = Ui=i 
(3.2) T:V —> S", T{x)=gi{x), x G V,. 

We need to find an open neigbourhood of the limit set where T is strictly expanding in the 
following sense : F defined on V is said to be (strictly) expanding on V with respect to the 
metric || || if there exists 9 > 1 {9 > 1) such that 

In the case when T is the positive part of a Schottky reflection group, we can suppose that 
up to a conformal identification dV^ is a great circle of the sphere 5”. For the metric on we 
can take the metric induced by its embedding in The inversion, ae, is the restriction to 

S'” of the symmetry on K."+^ with respect of the euclidean hyperplane countaining dVi, hence 
an isometry. Each inversion ai,i = 1 is expanding on the ball Hence, the map T 

is expanding on V = i U aeVi, strictly expanding on any open set precompact in V. 

However, the map T is not expanding on T) in general. To circumvent that we need to consider 
refinements, , defined by recurrence: 

=V 

VN = PI pAf-l ^ 

Each set is a disjoint union The collection of sets {'D^}i coincides with the collection 

^ rAr = { 7 er : | 7 | = iv}. 

Here | 7 | is the combinatorial length with respect to the system of generators {gi, ■ ■ ■ ^ 921 }■ In 
other words, 

= x~^^ ... x^^Vx^ if 7 = xi ... ccat , 
where 'Dxfj = T’i if xm = gi- The iterated map, , is defined on T)^, and 

=7, 

The map T is strictly expanding on T>^ for N big enough as explained in the following lemma 
(see Lemma 9.2 in Lalley for a similar result). 

Lemma 3.1. Let T a Schottky group and V, T defined as in ( |3.2| ). There exist an integer N >1, 
a metric || ||r, defined on T)^, and a real > 1 such that ||DT(w)||r > fi,w G . The metric 

can be taken analytic on . 


Proof. Let us recall that any Mobius transformation 7 of which does not fix the point at 
infinity, 00, has an isometric sphere Sj, and that 7“^ is strictly contracting on any compact 
subset of its exterior (the unbounded component of \ 5 '.y). The sphere Sj is centered at 7“^oo 
and if r^ is its radius, we have (see [p8|) 

r"^ \\x - v\\ 

(3-3) Il7a;-7y|| = T|- nil -n — rj: X,?/ G R'^ \ { 7 -^ 00 } 

||x — 7 -^oo|| lly — 7 oo|| 


Up to a conformal transformation, we can suppose that T is a subset of Conf(R”), with the 
point at infinity in its ordinary set. No non trivial element in T fixes 00 and, taking in ( 3 . 3 ) as 
x,y points in the upper half-plane 11 ”+^ (and the Poincare extension of 7 to R"+^), we deduce 
that the set of radii {r.y,7 G T} accumulates only at 0 . 

Eor 7 = xi... xtv, we have 


7 (^) ^ ^7 U T)x 2 ...X^ C . . . C T)x^_^xn U 
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hence there exists Nq such that the isometric sphere Sj is included in if I 7 I > Nq. For 
such a 7 , the interior of the isometric sphere S^-i is included in its exterior contains 

all the 'Dxi,Xi ^ , hence 7 is strictly contracting on 7 (^ 7 ) C 

the map 7 = is expanding on C We have just proved the existence of 770 > 1 

such that ||Z3T^“(?ii)|| > 'qQ,w e V^°, hence there exist constants C > 0,6* > 1 such that 
\\DTP{w)\\ > C0P,w €VP,p> 1. 

Taking an integer N such that C0^ > 1, we define on T>^ the metric (introduced by Mather 

@) 

iV-1 

\\V\\r=Y.\\DTP{w)V\\, VgT^V^, 

p—0 

which concludes the proof. □ 


Let §" be a Grauert tube of S", that is a complex n-manifold containing S" as a totally real 
submanifold (that is all we need). Let us then choose open neighbour hoods^ Di,i = l,...,djv 
of in By further shrinking, we can suppose that the open sets Di are mutually disjoint, 
and that the real analytic maps T and ||L>T||r extend holomorphically to D = with 

II.DT’llr 7 /? for some /3 > 1. The open sets Di can be chosen to be a union Di = of 

open sets, each one biholomorphic to the ball i?c"(0,1) in C". 

With this formalism in place we define the Ruelle transfer operator 


C{s)uiz) = ^ \\DT{w)\\ ®w(w), zGD, 


(3.4) 


Tw=z 


uGH'^{D), iL^(D) = {u holomorphic in Z) : JJ \u{z)\'^dm{z) < 00 } . 


The only difference from the standard definition lies in choosing spaces of holomorphic func¬ 
tions instead of Banach spaces. However we still obtain the analogue of a (special case of a) 
result of Ruelle I® and Fried [0 : 


Proposition 3.2. Suppose that C{s) : H^[D) H^{D) is defined by (3.4). Then for all s G C 

£(s) is a trace class operator and 

(3.5) I det(J — /l(s))| < exp((7|s|"'“'"^). 


Proof. The proof is based on estimates of the characteristic values, pi{C{s)). We will show that 
there exists C > 0 such that 

(3.6) /i^(/:(s)) 

To see how that is obtained and how it implies (^^) let us first recall some basic properties 
of characteristic values of a compact operator A : Hi ^ H 2 where L/j’s are Hilbert spaces. We 
define 

\\M = Zo(^) > ti-i{A) > ■ ■ ■ > pi{A) 0, 

to be the eigenvalues of [A*A)^ : Hi Hi, or equivalently of {AA*)i : H 2 i? 2 - The min-max 
principle shows that 

(3.7) fj,e{A) = min max \\Av\\h^ . 

codim V =■^ II u II jcjT^ =1 


^we drop the index N in the open sets for the purpose of notational simplicity 
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The following rough estimate will be enough for us here: suppose that {xj}'^Q is an orthonormal 
basis of Hi, then 


(3.8) 




J=t 


To see this we will use Vi = span {xj}"^^ in (0: for V G Ve we have, by the Cauchy-Schwartz 
inequality, and the obvious C inequality. 


||Ar 




'^{v,Xj)h,Axj 






y3=t 


from which (3.7) gives (3.8). 

We will also need some real results about characteristic values The first is the Weyl inequality 
(see |§, and also |^, Appendix A]). It says that if Hi = H 2 and \j{A) are the eigenvalues of A, 

|Ao(A)| > |Ai(A)| > ... > |A,(A)| ^0, 

then for any N, 

N N 

n(i+i^^(Ai)i)<n(i+i/^^(^)i)- 

e=o e=o 

In particular if the operator A is of trace class, that is if, '^ifJ-e{A) < 00 , then the determinant 


det(/+ A) — P(1 + A^(A)), 


e=o 


is well defined and 
(3.9) 


I det(/ + A)| < P(1 + l^i{Aj). 


e=o 


We also need to recall the following standard inequality about characteristic values (see |^): 
(3.10) /r^i+^2(A + i?) = ^^j(A) + ^^ 2 ( 7 ?) 

We finish the review, as we started, with an obvious equality: suppose that Aj : Hij —> H 2 j and 
we form 0 ^^^^ Aj : ^j^i Hij 0 /=i H 2 j, as usual, ®'j^i Aj{vi®-■ -©uj) = AiUi®- • •©Ajuj. 


Then 

(3.11) 


J 00 


E0“EE• 


f=0 




3 = 1 f=0 


With these preliminary facts taken care of, we see that (^) implies (|3.5|). In fact, (|3.9D shows 
that 

00 

,|n + l 


\ det{I - C{s))\ < P(l + e‘^l"l“^"/‘^) < 


i=e 

Hence it remains to establish (|3.6[ ). For that we will write 

dN 


H\D)^^H^{D,), 


2=1 


and introduce, for i^j = 1,..., the operator 


C,,{s) H\Di) ^ H\D^), 
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non zero only when T{Di) and Dj are not disjoint, where 

(3.12) Z&D,, = . 

From S) and a version of (ra we then have 

M^s)) < ^ max 2^[f/c](Aj (s)) • 

To estimate fj,k{Cij{s)), let us recall that Di was taken as an union of open sets Dik,k = 
1,..., (5i biholomorphic to Bc’^{0, 1): as fij{Dj) is relatively compact in Di, we can find p S (0,1) 
(independent of i,j = 1,.. .,dN) such that fij{Dj) C where '^ith C Afe 

the pullback of the ball i3c"(0, p) through the biholomorphism of Dik onto i?c"(0,1). The map 
Cij{s) is the composition 


HHDi) 










where R and ii(), are the natural restrictions, is the orthogonal projection on the space H^{D^) 
i mme rsed in by the natural restrictions and Cij{s) is defined by the same formula 

( 3.12 ) as Cij{s). The maps R and tt^ are bounded, while the norm of Cij{s) is bounded by 
The bounds on the singular values of given up to a bounded factor by the following 
lemma, give the bound 

for some C, which completes the proof of (|3.6D . □ 


Lemma 3.3. Let p G (0,1) and Rp : H^{Bc^{0, 1)) ^ H'^{Bc^{0, p)) induced by the restriction 
map of Bcn{Q, 1) to Bc^(0, p). Then, for any p G {p, 1) there exits a constant C such that 

MR^) < . 

with the standard basis (a:Q)Q,gN" of H^{Bc^{0, 1)): 


Proof. We use 
(3.13) 
for which we have 


,{z) =. 


/Bc"(0.1) 


\xa{z)f dm{z) = 1 , a G Ng 


A(^a)r = / \Xa{w)\‘^dm{w) = p 

J B^n (0,p) 


— ^2(1“!+") 


The number of a’s with |a| < m is approximately m” and hence by (3.8) we have 

j^n-lpk < 


Pi{Rp)<C Y. 

Ck>e/' 


□ 


The next proposition is a modification of standard zeta function arguments - see [ p^ and [ p^ 
for the discussion of the hyperbolic case. 

Proposition 3.4. Let C{s) he defined by (p.4|). Then, if Zy- is the zeta function (0) corre¬ 
sponding to the group T then 


Zr{s) = det(/ — /l(s)). 
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Proof. For s fixed and z G C 


h(z) det(/ — z£(s)) 

an entire function of order 0. For \z\ sufficiently small log(/ —z£(s)) 


is, in view of ( |3.6| ) and 
is well defined and we have 


(3.14) 


det(/ — zjC(s)) = exp ( — — tr(>C(5)) 

' n 


n—l 


The correspondence between the closed geodesic (or, equivalently, conjugacy classes of hyper¬ 
bolic elements) and the periodic orbits of T is particularly simple for Schottky groups and we 
recall it in the form given in (where it is given in a more complicated setting of co-compact 
groups): 

Closed geodesics on r\IHI, 7 of length /(y), and word length jyj are in one to one 
correspondence with periodic orbits {x,Tx,--- ,T'^~^x} such that ||T”(x)|| = 
exp^( 7 ), and n = lyj. For prime closed geodesics we have the same correspon¬ 
dence with primitive periodic orbits of T. 

It is not needed for us to recall the precise definition of the word length. Roughly speaking it is 
the number of generators of F needed to write down 7 . 

To evaluate tr £(s)™ we write 


trC{sr= tr(A 




where in the notation of ( p.l 2 | ) we have 

Chn O O A^_ii„M(z) = ||£>(/ni2 ° ° fim-limiz)) , 

o • • • o 


/u*2 o • • • o fi.^_.^i.^ : dq 




The trace of this operator is non-zero only if o ■ ■ ■ o has a fixed point in Di^. Since 

this transformation corresponds to an element of F that fixed point is unique. Let us call this 
element 7 “^. Since it corresponds to a given periodic point, x, of T", (corresponding to a fixed 
point of /iji 2 o • • • o 7 is determined uniquely by x and n: 

7 = 7 (x, n ), T”x = X . 

By conjugation and a choice of coordinates z it can be put into the form 


j{z) = 






and the trace can be evaluated on the Hilbert space 7L^(ilc"(0,1)). Using the basis (3.13) we 
can write the kernel of £ 7^2 o • • • o 


£722 o * • * o £j^_,,7 


Xz,w) = | 7 '( 0 )| ® ^ c „(7 ^(z)Y 


= E 


aGNJ 

(. g-(s+|a|)t(7)-i(S(7),«>^a^« 


aGN" 


The evaluation of the trace is now clear. 
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Returning to (3.14), we obtain for Re s sufficiently large (using { 7 }’s to denote the conjugacy 
classes of primitive of elements of F), 


det(/-z/:(s))=exp ^ ^ 


-{s+|o;|)^(7(rc,n))-2(0(7(a;,n)),a) 


n 

n—l T'^x—xocGNJI 


/ OO CXD 

= exp -EEEV^' 

n—l {-y} k—1 

\ h\=-n' 


-k({s+\a\)e{'l)-i(9('f),a)) 


=nn ( 1 _ ^l7lg-*{®(7).«>g-(s+|a|)^(7) 

{7}«GNJ 

which in view of o proves the proposition once we put z = 1. 


□ 


Remark. The proof above is inspired by the work on the distribution of resonances in Euclidean 
scattering - see |^, Proposition 2]. The Fredholm determinant method and the use of Weyl 
inequalities in the study of resonances were introduced by Melrose [|^ and developed further by 
many authors - see and references given there. That was done at about the same time 

as David Fried (across the Charles River from Melrose) was applying the Grothendieck-Fredholm 
theory to multidimensional zeta-functions [Q. In both situation the enemy is the exponential 
growth for complex energies s, which is eliminated thanks to analyticity properties of the kernel 
of the operator. 

Finally, we remark that in view of the lower bounds on the number of zeros of obtained in 
in dimension two, and in ||^ in general, we see from Proposition 3.4 that the upper bound 
3.5|) is optimal for any 7 . 

4. Quasi-self-similarity of A(r) and Markov partitions for r\IH”+^ 

In this section we will review the results on convex co-compact Schottky groups coming es¬ 
sentially from and [p^. The geometric point of view presented here was explained to the 
authors by Curt McMullenjj. 

We start with a more general definition of convex co-compact subgroups of Isom(IHI"’+^). A 
discrete subgroup is called convex co-compact if 

(4.1) r\C'(A(r)) is compact, (7(A(r)) convex hull(A(r)) 

Here, the convex hull is meant in the sense of the hyperbolic metric on 11"+^: A(r) C 

and 7 acts on it in the usual way. In particular this implies that r\C'(A(r)) has a compact 

fundamental domain in 11 "+^. 

The first result gives a quasi-self-similarity for arbitrary convex co-compact groups: 

Proposition 4.1. Suppose that T C Isom(IHI"+^) is convex co-compact in the sense of (i3)- 
Then there exist c > 0 and rp > 0 such that for any S A(r) and r < ro there exists a map 
g : R§n (xq , r) — » S" with the properties 

g(A(r) n Rsn(a;o,r)) c A(r) 


(4.2) 


cr ^dsr.{x,y) <d§r.{g{x),g{y)) <c V ^d§r.{x,y), a;,y G Rgn( xq, r). 


^who also generously provided us with Figures UJa, and previously with the essential part of Figure I 
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Proof. We proceed following the argument in |^, Sect. 3]. Let us fix zq £ C'(A(r)). If L is the 
geodesic ray through zq and xq, then 

3C>0yzGL3"fGT d{'y~^zo, z) < C. 


This follows from the compactness of r\C'(A(r)): for any point on the ray, z, there exists an 
element of the orbit of zq within a finite distance from z. We can now choose z = z(r) on the 
ray L so that d(z, zq) = log(l/r), and then 7 such that d( 7 “^zo, zq) = log(l/r) + 0(1). 

If x.y is the end point of the geodesic ray through zq and 7 “^zo, then for a fixed Ci, the 
ball Cir) covers B%n.[xQ,r). The action of 7 on Sg"( 2 ^ 7 , Oir) satisfies (^ ): A(r) is 

T-invariant, and the other property follows by putting 7 into the normal form (2.1). Since zq 
was fixed and we have no dependence on 7 , the proof is completed. □ 


Strictly speaking we will not use the quasi-self-similarity explicitely. It is present implicitely 
in the proof of Proposition 4A below where again the com pactness of r\C'A(r) is crucial. We 
remark that for n = 1, that is for surfaces, Proposition AA shows that, in the upper half plane 
model chosen so that A (s R, 


N(h) 

^(r) + [—/i, ^] = [o-jW,bj{h)], bj{h) < aj+i{h), \bj{h) - aj{h)\ < Kh , 
i=i 

with K independent of h. That gives a simple proof of our main result - see Section 

For the general we recall the definition of a Markov partition in a simple form |Q applicable 
here. Thus let /iA be the d-Hausdorff measure restricted to A(r). Let T and V be defined as in 

(ED- 

A Markov partition for the Schottky group P is given by a collection, V, of subsets oiV cEP, 
satisfying 

yP,Q€V tiA{T{P) n g) ^ 0 => T{P) D Q 
T is a homeomorphism on a neighbourhood of P n {T\p)~^{Q) if pla{T{P) CiQ) 0 

yPeV HAiP)>0 

VP,geP,p^Q /iA(pnQ) = o 
PAiT{P)) = Y. 1^a{Q) ■ 

QG-P 

fi^iT{P)nQ)^0 

A refinement of P is a new Markov partition given by 

n{V) = {{T\p)-\Q) P,Q€V, MA(P(P)nQ)^0}, 

see Fig.||. The refinements can of course be iterated giving PMifP). 

The Markov partition is said to be expanding if there exists a metric || || on P and a real 
9 > 1 such that ||P>r(u)^|| > 6>||^||,ii TyV. An example of a Markov partition for the 

Schottky group map is given by 

(4.4) ro = {Vi : l<i<2£}. 

If N is big enough its refinement, considered in Section 4, is expanding. 

The result we need to estimate the zeta function is given in the following 


Proposition 4.2. Let T be a eonvex co-eompaet Schottky group as defined in Seetion 2. Then 
there exist positive constants hg, Cq, and Ci, such that for 0 < h < ho there exists a Markov 
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Figure 5. A Markov partition and its refinements for the reflection group Fg 
of Fig. 1 (considered here as a Kleinian conformal group on 9IHI^ ~ U oo). 


partition, V^h), /or F with the following properties: 


(4.5) 


y P ev{h)3xp , PcSsn(xp,Co/i) 
yp,QGrih) A(F) nr(p) nQ 0 ^ dsni{T\p)-^{Q),dP) > h/Ci. 


In fact, we can make each element of'P{h) a union of a fixed number of balls of radii comparable 
to h. 


Proof. Rather than use the successive refinements of the Markov partition (which would 
not work) we will apply the following geometric observation: the projection of any geodesic in 
C'(A(F)) to F\IHI"+^ intersects the projections of C{d'Di) at uniformly bounded intervals, say by 
do > 0. This follows from the bound on the diameter of the compact fundamental domain of 

F\C(A(r)). 

For e < h < e we choose 'P(/i) in three steps. Let zq £ be a fixed point, and Vq 
the partition . We then put 

Vk(h) = {P &'R’^{Vo) : m < dw^+i{C{P),zo) <m +do] 

V{h) =' U ^ ■■ P(fQ& Vk-i{h)} . 

k>l 

The choice of do shows that the elements of P{h) cover A(F). From d{C{P),zo)) 
we conclude that P C i?sn(a;p. Cod), Co uniformly bounded, h ^ e“™. 

To obtain the strong form of the Markov partition property given in the second 
we modify V{h) as follows: 

(4.6) :P(d) = {{a: : d(x, A(F) n P) < d/Ci} : P £ P(d)} . 

We still have the property that 

V P, Q £ V{h) P^Q 


= m + C(l), 


part of (4.5) 


PnQnA(F) = 0. 
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If V(h) 5 Q C Dj, where Dj’s are as (3.12) then fij is a contraction of Q into a V(h) 9 Q' C Di. 


In fact, fij is a contraction of Dj into Di and it preserves A(r). That shows the second part of 


(4.5). 


Depending on the contraction constants of /ij’s (a finite set) we can replace elements of P(h) 
with unions of a hxed number of balls (the number depending only on the constants Cq and Ci 
and the dimension) so that all the properties still hold. □ 

5. Estimates in terms of the dimension of A(r). 

In Proposition |3.4| we used the Markov partition (4.4). It is clear from the proof that a Markov 
partition for which (4.5) holds (for some fixed h) can be used instead as long as the contraction 
property is valid in some complex neighbourhood of the components of the partition. For the 
proof of Theorem stated in Section 1 we will modify Df s in the definition of C{s) in a way 
dependent on the size of s: we will use Proposition ( p~^ ) with h = l/|s|. The self-similarity 
structure of A(r) will show that we can choose Dj = Dj{h) to be a union of disjoint 

balls of radii ~ h. A modification of the argument used in the proof of Proposition |3.2| will then 
give ([]|). 

We start with the following lemma which is a more precise version of the argument already 


used in the proof of Proposition 3.2 


Lemma 5.1. Suppose that Dj C C", j = 1,2, are open sets, and Di = Bc^{zk,rk)- Let g 
be a holomorphic mapping defined on a neighbourhood, Di of Di with values in D 2 , satisfying 

dc" (g(Di), 5 D 2 ) > I/Cq > 0 . 

If 

A : If^(D 2 ) Au(z) =\(g(z)), 

then for some Ci depending only on r^ ’s, K, dc’*{Di,dD 2 ), and min^^ ||I?5||c"^C"> we have 

t^i{A) < , 

where pf,(A)’s are the characteristic values of A. 

Proof. We define a new Hilbert space 

K 

Bk = BcAzk,rk), 

k^l 

and a natural operator 

J : H^iDi) n, {Ju)k = u\B, . 

We easily check that J*J : id^(Di) ^ H^{Di) is invertible with constants depending only on K. 
Hence 


ge{A) = n,{{rj)-^rjA) < ||(JV)-i|| \\r\\ni{JA). 


We then notice that 


gki{JA) < k max p.t{Ak ), 
l<k<K 


where 


Ak : H‘^{D 2 ) —> H‘^{Bk), Aku{z) = u{gk{z)), gk=9\Bk 


To estimate the characteristic values of Ak we observe that we can extend gk to a larger ball, 
Bk (contained in Hi) and such that the image of its closure still lies in H 2 (since we know that 
min^^ ||Z)g||c"^C" is strictly less than 1). That gives us the operators Rk : H‘^{Bk) H‘^{Bk), 
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RkU = and Ak defined as Ak but with Bk replaced by Bk- We now have Ak = RkAk and 

consequently, 

(J-eiAk) < \\Ak\\fJ.e{Rk) ■ 

Lemma 3.3 gives ^J,e{Rk) < (72 exp(—completing the proof. □ 


Proof of Theorem. As outlined in the beginning of the section we put h = l/|s|, where |s| is 
large but | Re s| is uniformly bounded. In Proposition |4^ each element of the Markov partition is 
given as union of (a fixed number of) balls. We can complexify this set by taking a corresponding 
union of balls in C" and all the properties hold for the analytic continuations of fij ’s defined in 

The now classical results of Patterson and Sullivan on the dimension of the limit set show 
that the total number of the balls is 0{h~^): what we are using here is the fact that the Hausdorff 
measure of A(r) is finite. 


We can now apply the same procedure as in the proof of Proposition 3.2 using Lemma 5.1. 
What we have gained is a bound on the weight: since | Re s| <C and is real on the real 

\[fij{z)Y\ < Cexpdsll arg/b(z)|) < (7exp((7i|s|| Itoz|) < C 2 , z€ Dj{h). 

We write C{s) as a sum of fixed number of operators Cij[s) each of which is a direct sum of 
0{h~^) operators. The balls and contractions are uniform after rescaling by h and hence the 
characteristic values of each of these operators satisfy the bound fie < C"f\ 0 < 7 < 1. Using 


(3.9) and (3.11) we obtain the bound 

log I det(/ - £( 5 ))I < CP{h) = 0{h -^), 

and this is ©. □ 

Proof of Corollary 1. The definition of Zt{s) dO) shows that for Res > (7i we have |^r(s)| > 
1/2. The Jensen formula then shows that the left hand side of ( |1.3D is bounded by 

y^{mr(s) : \s — ir — CA < € 2 } max log |^r(s)| + (£4 , 

Z-/ l^l^r+Cs 

Res|<Co 

and (|l.3| ) follows from (^). □ 


6. SCHOTTKY MANIFOLDS AND RESONANCES 

We recall that s complete Riemannian manifold of constant curvature —1 is said to be Schottky 
if its fundamental group is Schottky. In low dimensions Schottky manifolds can be described 
geometrically. 

Proposition 6.1. Any convex co-compact hyperbolic surface is Schottky. 

This result is proved by Button Q and for the reader’s convenience we sketch the proof. 

Proof. Any convex co-compact surface is topologically described by two integers {g, f) : its 
numbers g of holes and / of funnels, with the conditions <? > 0, / > 1 and / > 3 if 5 = 0. For 
any such pair {g, /), there does exist a Schottky surface of this type and we choose for each type 
(g, /) such a surface Mgj. The projection onto Mgj of the boundary of the Schottky domain is 
a collection £ 1 ,..., £^, of mutually disjoint geodesic lines. 

Let M be any hyperbolic convex co-compact surface . The surface M is homeomorphic to 
some Mgj. Pushing back on M the geodesic lines = Mgj and cutting M along 

these curves, we obtain in the hyperbolic plane a domain whose boundary is the union of paired 
mutually disjoint curves Ci,Ce+i,i = 1,... ,f, each one with a pair of points at infinity. These 
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points pair determine intervals, which are mutually disjoint (the curves Cj,j = don’t 

intersect). The intervals are paired with an hyperbolic transformation, so give a Schottky group, 
which coincide with the fundamental group of the surface M. □ 


Proof of Corollary 2. For a Schottky manifold, the fundamental groups is Schottky, and hence, 
M = r\]HI"+^,r C Isom’*'(IHI"’+^). We then introduce its zeta function Zm as the zeta function 
Zr of the group F. Following Patterson and Perry we introduce the spectral sets Vm and 
Sm defined by the Laplace-Beltrami operator Am on M: 

Vm = {s : Re s > n/2, s(n — s) is a Lf eigenvalue of Am} 

Sm = {s : Re s < n/2, s is a singularity of the scattering matrix Sm } 


Moreover, each complex s in Vm has a multiplicity denoted by rriMis), each s in Sm a pole 
multiplicity denoted by In the case of surfaces, the divisor of the Selberg zeta function 

Zm is given by the following formula: 


OO 

-Xm + l)[-k] + tum 



+ tom(s)[s]+ ^ m^(s)[s], 
s^Sm 


where xm is the Euler-characterictic of M, see Theorem 1.2] The zeta function, Zm, is 
entire and in any half-plane {Res > —Cq}, and the formula above shows that the bounds on the 
number of its zeros provide bounds on the number of resonances. The dimension of the limit set, 
S depends only on P and, as shown in , it gives the Hausdoroff dimension of the recurrent set 
for the geodesic flow on T*M by the formula 2(1 + 6). □ 

For a convex co-compact hyperbolic manifold M, Patterson and Perry give a formula for the 
divisor of the zeta function Zm in any (even) dimension, but it does not imply (in the non- 
Schottky case) that the zeta function is entire. In the case of Schottky groups, the zeta function 
Zm is entire, as it was shown in the Proposition ^.4 Hence we concluded that Corollary 2 holds 
also for Schottky manifold. 


We conclude with some remarks about Kleinian groups in dimension n + 1 


3-manifold are geometrically described by Maskit 114 


3. Schottky 


Proposition 6.2. A hyperbolic convex co-compact, non compact 3-manifold is Schottky if and 
only if its fundamental group is a free group of finite type. 


While non compact surfaces of finite geometric type have always a free fundamental group, 
that is not true for 3-manifold. For instance, if F is a co-compact surface group, the 3-manifold 
H^/F is convex co-compact with a non-free fundamental group. Quasi-fuchsian groups (that is, 
deformation of such a F in Isom(IHI^)) give similar examples. 

Finally, we note that the bound on the number of zeros of Zr established here for Schottky 
groups is valid for any group F, for which a expanding Markov partition can be built. Anderson 
and Rocha construct such Markov partition for any function group. This class of groups 
does not exhaust all convex co-compact groups (the complement in the 3-sphere of a regular 
neighbourhood of a graph is not in this class) and it is not known if all convex co-compact 
Kleinian groups admit an expanding Markov partition. 


7. Numerical Results. 

7.1. Discussion. In this section, we present numerical results on the distribution of zeros of a 
closely related dynamical zeta function Z(s) for a conformal dynamical system. We consider the 
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9 

dim 

10 ° 

0.11600945 

20 ° 

0.15118368 

30° 

0.18398306 

0 

0 

0.21776581 


Table 1. Dimensions of the limit set for relevant values of 9. 


simple case of groups Fg generated by reflections in three symmetrically placed circles perpen¬ 
dicular to the unit circle (considered as the boundary of the Poincare disc), and intersecting it 
at angles 9 - see Fig.^ where 9 = 110°. 

Numerical computations of the zeta function in that case have been already performed by 
McMullen and Jenkinson-Pollicott |^. Their goal was to find an efficient way of computing 
the dimension of limit sets. Table ^ gives the (approximate dimensions of the limit sets for the 
relevant angles, calculated as the largest real zero of Z{s) Q using Newton’s method. 

Figures |7[^ show 

^ log(|{s g [xo,xi] X [yo,y] : Z{s) = 0}|) _ 

^ ^ log(2/) 

as a function of y. In each plot, the value of Xq is varied to test the dependence of the distribution 
on the region in which we count: The blue line corresponds to xq = —0.2, the red line xg = —0.1, 
and the black line xg = -kO.l. The data show that most of the zeros lie in the left half plane. 
Based on the theorems proved in earlier sections, we expect the curves to be bounded above by 
the dimension (the thick blue line) asymptotically. This is not the case, except for the black line, 
which represents zeros with Re(s) > xg = -1-0.1. 

Similarly, Figures show '°s(']°s(i|-^|WI)) function of |s|, for a large number of points 

in the rectangle [—0.2,1.0] x [0,10^]. In this case, we also expect the curves to be asymptotically 
bounded by the dimension. This is also not the case. The only reasonable explanation, barring 
errors in the numerical calculations, is that the asymptotic upper bound is accurate only for very 
large values of Ito(s), and we were not able to calculate Z(s) reliably for such values. These 
results also show that Z(s) has plenty of zeros in regions of interest. 


7.2. Implementation notes. To count the number of zeros of Z(s) in a given region fi in the 
complex plane, we rely on the Argument Principle: 

(7.2) \{,^n:Zis)=0}\ = :^ [ ^ds. 

2711 Jon Z{s) 

To evaluate Z{s), our main technical tool comes from Jenkinson and Pollicott Q, though we 
note that the essential ideas were used in Eckhardt, et. al. ||] and date back to Ruelle [ p^ . 

First, some notation: Let us denote symbolic sequences of length jyj = n by 7 . That is, 
7 = ( 7 ( 0 ), 7 ( 1 ),..., 7 (n)), 7 (fc) e {0,1,2}, and 7 ( 0 ) = 7 (n). Such sequences represent periodic 
orbits of length n. To each such sequence 7 we associate a composition of reflections = 
o o : F* 7 (o) ^ D^(gy As 4>-~f is a contraction of D.^(o) into itself, it has a unique fixed 
point z-y. 

It is shown in Jenkinson and Pollicott that Z(s) = limM —>00 Zm{s), where 


Zm{s) = 1 + ^ 


E n- E 

r! rife 

N=lr=l neP{N,r)k=l |7|=nfc 




(7.3) 









THE SELBERG ZETA FUNCTION FOR CONVEX CO-COMPACT SCHOTTKY GROUPS 


19 



Figure 6 . This plot shows \Z{s)\ in the square [—0.2,0.3] x [—2,2], for 9 — 30°. 



Imfsl 

Figure 7. This plot shows i°g(l{^s[^o,a:i]^x|yo,y]:Z(s)=o}|) _ a,s a function of y, 
for different values of xq\ The thin blue line is for xq = —0.2, the red line for 
Xo = —0.1, and the black line for Xq = +0.1. Note that the value of xi is not 
very important because Z{s) decays very rapidly for large Re(s). Thus, we set 
xi = 10 throughout. The value of yo is fixed at —0.1, to avoid integrating over 
any zeros. The thick horizontal line indicates the dimension of the corresponding 
limit set. In this plot, 9 = 10°. 
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Figure 8. This plot shows oil) — 1 as a function of y, 

for different values of xo- The thin blue line is for xo = —0.2, the red line for 
Xo = —0.1, and the black line for xo = +0.1. Note that the value of xi is not 
very important because Z{s) decays very rapidly for large Re(s). Thus, we set 
xi = 10 throughout. The value of yo is fixed at —0.1, to avoid integrating over 
any zeros. The thick horizontal line indicates the dimension of the corresponding 
limit set. In this plot, 9 = 20°. 


where P(N, r) is the set of all r-tuples of positive integers (ni,..., Ur) such that ni + ... + nr = N. 
The series (in N) converges absolutely in {s : Re(s) > —a} for some positive a. 

Equation ( |7.3| ) lets us evaluate Z{s) for reasonable values of s in a straightforward manner. 
However, we found two simple but useful observations during the course of this calculation: 

(1) Define 


(7.4) 


a„(s) = 


1 


and 


(7.5) 


BN,r{s) 


n^P{N,r) k—1 


(7.6) 


Then the recursion relation 


BN,r{s) = -- ^2 BN-n,r-l{s) ' ara(s), 


with initial conditions B^ i(s) = aAr(s), provides an efficient way to evaluate the sum in 
A similar relation can be derived for Z'(.s) by differentiation. 

(2) Recall that the maps 4>.y are compositions of linear fractional transformations. Identifying 
such maps with matrices in G'L(2,R) in the usual way, we can compute the numbers 
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lm(s) 


Figure 9. This plot shows oil) — 1 as a function of y, 

for different values of xo- The thin blue line is for xo = —0.2, the red line for 
Xo = —0.1, and the black line for xo = +0.1. Note that the value of xi is not 
very important because Z{s) decays very rapidly for large Re(s). Thus, we set 
xi = 10 throughout. The value of yo is fixed at —0.1, to avoid integrating over 
any zeros. The thick horizontal line indicates the dimension of the corresponding 
limit set. In this plot, 9 = 40°. 


via matrix multiplications. However, such matrix multiplications can become 
numerically unstable for larger values of jyj. 

An alternative involves the observation that the matrices cor¬ 

responding to the maps (f>j have distinct nonzero real eigenvalues. Let us denote these 
eigenvalues by A+ and A_ so that |A+| > |A_|. Then a simple calculation shows that 
(jj'-yiz'y) = A_/A+. This becomes simply (—1)I'>'I/A^ if we normalize the determinants of 
the generators Aq, Ai, and A 2 . The larger eigenvalue A+ can be easily computed using 
a naive power method: 

(a) Choose a random vq. 

(b) For each A: > 0, set Vk+i = A.yVk/\\A.yVk\\ and A^^ = {AryVk,Vk). 

(c) Iterate until the sequence (A^^) converges, up to some prespecified error tolerance. 
The resulting algorithm is slightly less efficient than direct matrix multiplication, but it 
is much less susceptible to the effects of round-off error. 

Note that it is certainly possible, even desirable, to apply to this problem modern 
linear algebraic techniques, such as those implemented in ARPACK 0. But, we found 
that the power method suffices in these calculations. 


These two simple observations lets us calculate the values of Z{s) for a wide range of values in 
an efficient manner. When combined with adaptive gaussian quadrature. Equation (7.3) allows 
us to evaluate the relevant contour integrals. 
























22 


L. GUILLOPE, K.K. LIN, AND M. ZWORSKI 



Figure 10. This plot shows log(log(|.Z(s)|))/log(|s|) for a large number of 
points in the rectangle [—0.2,1] x [0,10^]. The horizontal line indicates dimen¬ 
sion. Here, 0 = 10°. 



Figure 11. This plot shows log(log(|.Z(s)|))/log(|s|) for a large number of 
points in the rectangle [—0.2,1] x [0,10^]. The horizontal line indicates dimen¬ 
sion. Here, 9 = 20°. 
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Figure 12. This plot shows log(log(|.Z(s)|))/log(|s|) for a large number of 
points in the rectangle [—0.2,1] x [0,10^]. The horizontal line indicates dimen¬ 
sion. Here, 6 = 40°. 



Isl 

Figure 13. Logarithmic plot (base 10) of the modihed relative error 
along the line Re(s) = -0.2, , where Rn{s) = Z'j^{s)/Zn{s). 
The blue curve is 0 = 10°, the red curve 9 = 20°, the green curve 9 = 30°, and 
the black curve 9 = 40°. 
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Figure 14. Logarithmic plot (base 10) of the modified relative error 
along the line Re(s) = -0.1„ where Rn{s) = Z'j^{s)/Zn{s). 
The blue curve is 0 = 10°, the red curve 0 = 20°, the green curve 9 = 30°, and 
the black curve 6 = 40°. 



Isl 

Figure 15. Logarithmic plot (base 10) of the modified relative error 
along the line Re(s) = +0.1, where Rn{s) = Z'j^{s)/Zn{s). 
The blue curve is 9 = 10°, the red curve 9 = 20°, the green curve 9 = 30°, and 
the black curve 9 = 40°. 
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Remarks. 


(1) To calculate the Selberg zeta function ^ 2 ( 5 ) for closed geodesics on the quotient space 
r\]HI^, we simply sum over periodic orbits of even length, and additionally use 02 (n, s) = 
2a(n, s) instead of a(ri, s) in the recursion relations above. This counts the number of 
equivalence classes of orbits correctly. 

(2) The work of Pollicott and Rocha revolves around a closely-related trace formula: 


(7.7) 


Z{s) 


00 N r 

i + EB-i)’' E n 

N =1 r —1 {[7i],...,[7r]}ePr(A',r) fc=l 


|1 



2 


where Pr(iV,r) = {{[ 71 ], [ 7 r]} : | 7 i| + ... + | 7 r| = primitive}, and [ 7 ] is the 

equivalence class of 7 under shifts. The primary difference between (7.3) and ( 0 ) is that 
the latter sums over sets of equivalence classes of primitive periodic orbits (equivalent 
up to shifts), whereas the former sums over all periodic or bits. While it is possible to 
enumerate primitive periodic orbits efficiently,! Equation (O) still provides a better 
numerical algorithm, as it is easier to implement and results in faster and more stable 
code. 


7.3. Error analysis. Figures |13|-|15| show the logarithms (base 10) of the modified relative errors 


/y ox |Ri2(s) - Ri3{s)\ 

l+|Rl2(s)| + |Rl3(s)r 

on the lines xq -I-j[0,10^], for xq S {—0.2,—0.1,0.1} and where Rn{s) = Z'jy{s)/Zpf{s). This 
formula interpolates between the absolute and the relative errors, and measures the convergence 
of the integrand in (([7. 2 ])). These results lend some weight to the reliability (i.e. convergence) of 
the values of Z'j^{s)/Zn{s) used in the calculations above. 
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